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Abstract 



An application of a self-consistent version of RPA to quantum field theory 
with broken symmetry is presented. Although our approach can be applied to 
, any bosonic field theory, we specifically study the theory in 1+1 dimensions. 

CS| ' We show that standard RPA approach leads to an instability which can be 

removed when going to a superior version, i.e. the renormalized RPA. We 
' present a method based on the so-called charging formula of the many electron 

, problem to calculate the correlation energy and the RPA effective potential. 

1 Introduction 

Oh! 

^ ■ One central aim of the relativistic heavy ion program is to produce highly excited hot 

and dense matter possibly constituting a quark-gluon plasma. In such a phase quarks 
and gluons should be liberated and chiral symmetry should be realized in its Wigner 
^ ■ form. A central theoretical question is thus to have a correct description of the 

broken vacuum and of the progressive restoration of chiral symmetry with increasing 
temperature and/or baryonic density. This problem is highly non-perturbative in 
nature and usual perturbative loop expansion technique are certainly not sufficient 
by construction. These features provide at least one important motivation to develop 
tractable non perturbative methods to be applied in the context of (effective) field 
theories with broken symmetries. 

One hope is to try to adapt to quantum field problems very well controlled non 
perturbative methods from the nuclear many-body problem possibly of variational 
nature. A very popular method, known as the Gaussian approximation for in- 
teracting bosons, exactly corresponds to the Hartree-Fock-Bogoliubov mean field 
approximation which constitutes the basic building block of the nuclear many-body 
theory i.e. for fermions, see e.g. [1]. This variational method has already been 
applied to theories with a global symmetry [2] as well as with a local symmetry 
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[3]. Nevertheless when this HFB approach is applied to a bosonic 0{N) model {i.e. 
to the linear sigma model) the pion appears with a finite mass or, in other words, 
Goldstone theorem is violated as discussed in [4]. As demonstrated in recent papers, 
RPA fluctuations are able to generate this soft mode, known as the spurious mode 
in the context of nuclear physics [5, 6]. The RPA approach thus appears as a very 
promising technique to treat non-perturbative problems in the context of quantum 
field theory. However in its simplest form this method is not variational and its 
predictions are not always very well under control. In particular, it is well-known 
that standard RPA has the tendency to overestimate the attractive correlation en- 
ergy, at least in examples of nuclear physics. The purpose of this paper is therefore 
two-fold: first, a complete presentation of the RPA technique is given in the context 
of the simplest field theory, namely the A<y9^ theory with a specific application to 
the 1-1-1 dimensional case. The second goal is to develop the formalism of more 
elaborated versions of the RPA approach, namely the so-called renormalizcd RPA 
(r-RPA), which, to our knowledge, has never been done before. Let us mention 
that there exist previous attempts to apply many-body techniques for this specific 
problem. For instance, Hauser et al. use the cluster expansion of Green's functions 
[7]. Of course other methods such as lattice calculation have been employed, see 
e.g. [8, 9, 10]. Our aim is not really to compete with these numerical methods 
but actually to develop a tractable approach allowing direct physical interpretations 
in view of further applications in the context of chiral or gauge theories. In this 
preliminary work our aim is rather modest and limited to the presentation of the 
formalism and the discussion of the remaining problems to be solved such as the 
explicit covariance of the obtained results. We also present the method to calculate 
the RPA correlation energy in the context of a field theory within a Green's function 
formalism. One important point of this paper is to adapt the so called "charging 
formula" method (see e.g. [11]) for quantum field theory beyond the standard RPA 
scheme. Some numerical results are also obtained showing in particular how to cure 
the instability of the standard RPA in the broken symmetry region. 

2 The theory 

2.1 The Hamiltonian 

We consider the Lagrangian density : 

£ = i d'^cpix) dM^) - i cp^x) - A/(^) (1) 

where /Xq is a constant and the bare coupling constant 6 = A/6 is positive for reasons 
of stability. We decompose the scalar field (p{x) in a classical part or condensate s 
and a fluctuating piece (/>(x) : 

(p{x) = (t){x) + s, s= {(p{x)) . (2) 



The presence of the condensate s indicates a spontaneous breaking of the underly- 
ing (f —ip symmetry. Introducing the conjugate field n(a;), one obtains for the 
Hamiltonian (in d+l dimensions) : 

H= jd'^X {1^^^^ + A,4 + ^^^^ + ^,3^ 

+ <t>\x) + ^^ct>\x)Y (3) 

Putting the system in a large box of volume = L'', it is convenient to work in 
momentum space and to expand the fields according to : 

'^(^) = ^T. m ' n(x) = - n,KO • (4) 

The hermiticity of (/>(x) and Ii{x) imposes = (f)_^ and lit = — H.^- and canonical 
equal-time commutation relations translate into : 

[•/-g-'nt,] = . (5) 

The Hamiltonian can be rewritten as : 
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+ ^ X] ^1,2,3,4 ^1 02 <?^3 '?!'4 (6) 
1234 

where (^j is a short-hand notation for and the three-body and four-body inter- 
actions are given by : 

Vl,2,3 = %+g2+g3' ^1,2,3,4 = y %+g2+g3 + g4 • (7) 

The bare single particle energy appearing in the Hamiltonian is : 

Ol = t + nl + \sK (8) 

In the following we will call and i?4 the three-body and four-body interacting 
Hamiltonians. 



2.2 The Gaussian approximation 

In the Gaussian approximation the ground state is represented by a trial wave func- 
tion which is a functional of the field ip{x) : 

m<p)) = AT exp (-i I d^x d% {^{x) - s) h{x - y) {ip{y) - s)^ . (9) 

The optimal wave function is obtained by minimizing {'^(ip)\H\'^{}p)) with respect to 
h{x — y). The resulting function, which still depends on the condensate ,s, defines the 
Gaussian effective potential. The various minima in s correspond to the possible 
phases of the system. Working in momentum space, one introduces the Fourier 
transform of h(x — y) according to : 

h{x -y) = j d'^Q. £q ■ (10) 

It is easy to show that the trial ground state is the vacuum of the canonical destruc- 
tion operator such that : 




The single particle excitation of this vacuum have energies Eq which differ from 
the bare energies Oq. In other words we have rotated the original bare basis with 
single particle energies Oq into a basis associated to the e^'s through a Hartree-Fock- 
Bogoliubov (HFB) transformation. In this HFB ground state the energy density is 
easily calculated by using the Wick theorem : 

«.,.) = i.§.^ + ^.' + ^Z(| + |)n<*V (12) 

where the scalar density {(fP') is given by : 

(^^>4s<#4) 4^4-/(04- (13) 

Minimization with respect to Sq gives the HFB quasi-particle mass m : 

m' = el-t=l4+\ is' + (02)) . (14) 

It has been demonstrated in d = 3 spatial dimensions that the above gap equation 
can be rendered finite by appropriate mass and coupling constant renormalizations 



[12]. Here we concentrate on to the d = 1 case where the theory becomes super- 
renormahzablc and only requires a mass renormahzation. We ehminate /xq in favor 
of the renormahzed mass /j, according to : 
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where A is a ultraviolet cutoff. The gap equation becomes 



m' = f^'+^- {s' + {A4>'),) (16) 



with 



which is independent of the cutoff A. Re-injecting the solution for the mass m 

into the expression of the energy density one gets the effective potential which is 
also finite, once the energy of the perturbative vacuum of particles with mass /x is 
removed. The result is : 

Snis) 1 9 A \ ( m? 'rin?\ 3p / , , 

where 

p = 6/2V 

is the dimensionless coupling constant. At low p the effective potential has a min- 
imum at s = : the unbroken phase is thus stable. Increasing p, a new minimum 
develops at finite s corresponding to a meta-stable broken phase. At a certain crit- 
ical pc = 2.57 the deformed phase becomes stable and the system undergoes a first 
order phase transition (see figure 1) in agreement with [7]. The approach can be 
straightforwardly extended to finite temperature using the so-called statistical vari- 
ational principle. One obtains for the finite temperature effective potential {i.e. the 
grand potential) : 

no(s,r) 
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In-^ -3p{Acl>X(.Tm 



with 



where n{x) = l/{exp{x) — 1) is the Bose-Enstein distribution and the gap equation 
(16) is unchanged once {A(f)'^)fj_ is replaced by its finite temperature expression given 
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Figure 1: Left panel : Gaussian effective potential for various values of the dimen- 
sionless coupling constant p. Right panel : value of the condensate s minimizing the 
Gaussian effective potential as a function of p. 



just above. Some numerical results are shown on figure 2. Choosing p > Pc the 
vacuum is in a broken phase. At a certain critical temperature Tc one gets a first 
order transition towards the symmetry restored phase. 



3 The RPA approach 



3.1 The equation of motion method (EOM) 

The aim of the RPA method is to describe the excitation spectrum of a Hamiltonian 
H. The excited states and the ground state |0) are defined by the conditions : 



Ql\o), 



Qu\o) = . 



(21) 



Minimizing = {vlHlv) /{vlv) with respect to the operators Qi^, one gets the 
following set of equations [1, 13] : 



(0| 



SQ,,[H,Ql] \0) = n,{Q\[5Q,, Qi]\0) 



(22) 



where i^i, = Ei^ — Eq is the excitation energy. Equation (22) is supplemented by the 
conditions : 

(0|[i7,Q,]|0) =0 (23) 
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Figure 2: Left panel : Gaussian effective potential for p = 3.8 for various temper- 
atures T in units of /i. Right panel : value of the condensate s minimizing the 
Gaussian effective potential as a function of T (in units of jj) for various values of p. 



which will determine the mean- field basis. Eq.(23) is a natural complement to 
eq.(22), since both equations hold in the exact case [14]. In general in this EDM 
approach, the excitation operators are searched only within a limited domain. They 
are usually taken in the form : 



(24) 



where the Aa constitute a set of non-hermitian operators labeled by a set of quantum 
numbers a (for instance momentum or isospin state). The RPA equations (22) now 
become matrix equations which allow to determine the excitation energy and the X 
and Y amplitudes (summation over repeated indices is understood) : 



with the matrix elements given by the double commutators : 



E* A* 



Yb" 



^b 



and the norm matrix is : 



|0), Sa6 = -(0| 



Aa,\ii,A^^\ 



|0) 



{0|[A„,4]|0) (0|[A„,^_,]|0) 
-(0|[^„,^_,]|0) -(0|[^„,4]|0)^ 



(25) 



(26) 



(27) 



The previous conditions (23) translate into equations having the form of gap equa- 
tions: 

{0\[H,Aa]\0)=0. (28) 

In practice the quaUty of the RPA scheme depends on the approximation which are 
made. Firstly the larger is the space of Aa operators, the better is the accuracy of the 
method. Here, in practice we will limit ourselves to single and pair boson operators 
namely A\ = {fot, fttfe^, ^^g*}- Secondly, one crucial point is the calculation of 
the ground state matrix elements of the double commutators and of the matrix 
elements appearing in the gap equations. Calculating the matrix elements on the 
true RPA ground-state constitutes the self-consistent RPA (SCRPA). In practice 
this is a formidable task both formally and numerically, which only has been fully 
achieved in a very limited number of simple model cases. In the most common 
approximation (i.e. the standard RPA), the ground state matrix elements entering 
the RPA equations are calculated on the mean-field ground-state. In our case it 
coincides with the Gaussian or HFB ground state discussed above. There is an 
intermediate approximation scheme, called in nuclear physics renormalized RPA 
(r- RPA) [13, 15, 16, 17] which allows to incorporate part of the correlations in the 
ground state. In section 4, we will propose the first application of this method in 
the context of quantum field theory. 



3.2 The Dyson equation approach (DEA) 

We present an alternative but equivalent formulation of the RPA scheme based on 
a Green's function (GF) approach [18]. Although equivalent in principle to the 
equation of motion approach, this DEA method is more convenient for quantum 
field theory. We define the time-ordered Green's functions : 

Ga,bit,t') = -^{0\T[Aa{t),Al{t'))\0) 

GaMt,t') = (Aait) , A_,it')) \0) 

G_,,-,(t,t') = -i{0\T (AUt) , Alit')) \0) 
G_,,_ft(t,0 = -i{0\T[Al,{t),A_,{t'))\0) . (29) 

We introduce the energy representation of these GF according to : 

G{t, t') = I ^ e-' ^ G{E) (30) 

and the matrix : 

G JF)-{ ^afii^) Ga-b{E) \ . . 



The RPA equations have now the form of a set of coupled integral equations for the 
various GF : 



E Gam = + E ( B^^^ ) 

„ J \ '^—a—c a— c / 



(32) 



with : 



Aab = (0| 



|0), 



ah 



-(0| 



|0). 



(33) 



Notice that the double commutators are ordered differently from the ones obtained 
in the equation of motion method. In practice, it can been shown that they are 
always identical. 



4 The RPA and the renormalized RPA appUed to the 
(^"^ theory 

4.1 Solution of the RPA problem 

We will limit ourselves to the case where the A operators are only one-body and two- 
body operators. We introduce creation 6^ and destruction operators bp depending 
on the parameters according to : 

or equivalently : 




By construction they obey standard canonical commutation relation for boson cre- 
ation and destruction operators : 

&^'] = <^/3,/3', [bp, V] = • (36) 

The /3's represent the quantum numbers of the created boson state; here this is 
simply a momentum index and —j3 represents the opposite momentum. 

The excitation operators A\ will be the one-body operators and the two-body 

operators 6g b]y and ftjj 6_^' with /3 / /?' for the "particlc-holc" operators. The gap 
equation {[H, hp]) = will give the extrema in the condensate s of the vacuum energy. 
Since we will calculate the effective potential giving directly the true minimum (i.e. 
the stable phase) we will not use it here. The most interesting non trivial gap 



equation is ([i?, = which constrains the basis, that is the parameters. 

This gap equation gives ((...) stands for ground-state expectation value) : 

^ 1,2 

+ ft 51 (<?^/3'/'l4</'3)conn ^1,2,3,-/3 (37) 
1,2,3 

where the sufHx conn means connected operators which contain exclusively corre- 
lated expectation values : 

{ABCD)^onn = {ABCD) - {AB){CD) - {AC){BD) - {AD){BC) (38) 

and £/3 is given by : 

e} = f,+^^l + \s'+\{<l>') (39) 

with : 




The Eps can be seen as the generalized mean field single-particle energies. But 
at variance with the Gaussian case these energies depend on the correlated Self- 
Consistent scalar-density (^^). 

Standard RPA 

In the standard RPA, we omit the connected parts. We thus get from the gap 
equation (37) : 

(n^nt ) - 4 ^<^^^t ^ = (41) 

and take for {(j)i3(j)p) its value in the Gaussian HFB ground state. This implies (see 
eq.(35) : 

= e/3 • (42) 

The basis is fixed and coincides with the HFB basis. Similarly all the matrix elements 
appearing in the calculation of the double commutators and the norm matrix are 
simply obtained by using Wick theorem on the Gaussian ground state. 

Renormalized RPA 

In renormalized RPA (r-RPA), one still systematically omits correlated expectation 
values. Hence the gap equation (41) remains valid. However, the scalar density is 
not yet fixed. It has to be determined self-consistently. The problem of its evaluation 
(in a non relativistic many-body problem for Fermi systems it corresponds to the 



occupation numbers) is one of the subtleties of self-consistent RPA (SCRPA). We 
will come to this problem later on. The gap equation can be rewritten as : 

(4 - 4) (1 + 2 ^'J bp) = (4 + a}) {bl hip + hp h^p) . (43) 

The basis np is not yet totally fixed. However, we have checked that the result of the 
r-RPA calculation does not depend on its choice, provided the above gap equation 
is satisfied. Thus we can choose it as up = ep i.e. the generalized mean-field basis. 
This choice has the merit of significantly simplifying the lengthy calculation of the 
double commutators. In that case one has : 

{hlhlp) = {hph_p) = Q. (44) 

In SCRPA this last property has to be always satisfied because the 6^ hip operators 
are just linear combinations of the genuine RPA excitation operators Q\, whose 
expectation values on the true RPA ground state vanish by construction. 
The many-body operator expectation values are also calculated using Wick theo- 
rem but the resulting two-body operator matrix elements depend on the various 
occupation numbers. In other words all the ground-state matrix elements enter- 
ing the double-commutators and norm matrices are expressible in term of the fixed 
momentum densities Mp = {(ppcp^p). We display here some examples : 

(hi>-D) = 4^^" 

(6146364) = (6l6i) (6^62) (5l,3 ^2,4 + ^1,4 ^2,3) 

+(6l6Li) (636-3) ((5i,_2 (53,-4) . (45) 

We do not give the details of the calculation to obtain the solution of the RPA 
problem (32). Even if we limit ourselves to one- an two-body operators, it represents 
quite a lot of algebra. The results for the various GF are listed in the appendix. 
Here we give only the GF relative to the field operator (f)p. For the one-particle GF 
i.e. the (f) particle propagator, one obtains : 

G<t>p,<t,p,{E) = 6pp,G{E,P) 

G{E,P) = (^E^ -s%- S(£;,P))"' (46) 

where the mass operator has the following form : 



62^2 _ I{E,P) 

1 - \I{E,P) 



j:{E,P) = ^i{E,P) with I{E,P)= ^ (47) 



and the two-particle loop has the explicit expression : 

^ "ei + £2 

riA'i — C2A2 



mP) = 



dki dk'}, 
ei - £2 



eiNi + g2A/"2 

2£ie2 - (£1 + e2)Vir? 



(48) 



Notice that correlations are present in this expression through the densities Hi = 

{(picf)]). We will see in subsection 4.4 how to calculate these densities in r-RPA. 
For the lp-2p and 2p-2p GF we give the particular combinations which are directly 
relevant for the calculation of the effective potential : 



123 

^1,2-3,-46'^ 

1234 



dP ^ - 



, aa{E) = V 



dP 
(27r)° 



2I{E,P) + b 



P{E,P) 



1 



I{E,P) 



I{E,P) 



■HE,P), 



G{E,P) 



(49) 



The reader may check that these results have a very clear diagrammatic interpreta- 
tion (see figure 4). Prom these expressions, one can get the expectation value of the 
three-body Hamiltonian by taking the appropriate t' ^ t limit of the GF. 



(50) 



The correlated part of the four-body Hamiltonian is obtained with the same tech- 
nique but the uncorrclatcd GF has to be removed. As wc will see below the expres- 
sion of the correlated energy will involve the following quantity : 



1 



2AV 
1 



idE 

W) 

dp 



1234 



24 J (27r)^ 



idE 
I{E,P) 



2,-3,-4 (G<^3^^^^t^t(^) - 

,2 IHE,P) 



M , <f>l4'l 



I{E,P) 



1 



mp). 



G{E,P) 



(51) 



where the first order term {i.e. the one loop term I{E,P) of eq. (49)) has been 
removed. 



4.2 The Effective Potential 

The starting Hamiltonian eqs. (6, 7) contains a free Hamiltonian of bare particles 
with energy and an interacting Hamiltonian + H4. To simplify the writing 
we now replace the momentum labels q by integers i and omit the linear term in (j) 
which does not directly play a role in the formal manipulations : 

H = v(^\tils^ + ^s^^ + Y^^{uM + Olh4) +H3 + H4. (52) 

In the RPA approach, the contribution of the interacting part of the Hamiltonian 
systematically transforms, in the expressions of the various double commutators, 
the bare single particle energy into the generalized mean-field single particle energy 
(39) which is finite after mass renormalization in one spatial dimension : 

Ol^el = 0!+^- {ct?)n (53) 

where {^^)r = (1/^) ^1) is in principle the correlated vacuum density (see 
subsection 4.4) except in standard RPA where it is taken on the Gaussian ground 
state. 

This suggests to rewrite the Hamiltonian as : 

H = Ho + Hint (54) 

with : 

Hint = Hs + H4-^ ^ : <^t - y ^ (02)2 _ (55) 

{^'^)e is the expectation value taken on the ground state of the mean-field quasi- 
particles with energies £\ and : ... is the normal ordering with respect to this 
vacuum, namely : 

■.(I)i<p{:e=^i4>\- {(l>i<Pi)e ■ (56) 

We use again the notation {4>'^)r for the scalar density. In renormalized RPA this 
scalar density is in principle the correlated one. In standard RPA the quantity 
appearing in Hint is identified with the mean-field scalar density ((^^)e and e refers 
to the Gaussian HFB mean-field. In that case, the interacting Hamiltonian reduces 
to : 

{Hint ) standard RPA = H3 + : H4 :^ . (57) 
The Hq Hamiltonian is obtained as H — Hint '■ 

+l{<i>')Rj:--'p^4--e+vi{4>r,. (58) 



Again in the case of standard RPA the self-consistent scalar-density (</>^)_r is replaced 
in the above expression by the Gaussian HFB scalar density {(t>^)e- Hq can be 
rewritten as : 

i^o = i^o + E ^ ni nl 4 -.e) . (59) 

It has a form of a free Hamiltonian for quasi-particles with mass m (see eq. (39)) 
i.e. w? = £^ — (f'. In one spatial dimension, this mass is rendered finite by a simple 



i2 = £2 _ ^ 
mass renormalization 



2 2,^2,^ ij.2\ f^^ dq 
m = iJ, + -s + 



2 2 J-A 27r 2V?+7?, 

Eq is the generalized mean-field vacuum energy : 

-^0 1 2 2 , ^ 4 I /^2\2 

V 2^° 24 8 



+ ^ 

1 



2^0' +24" + 2F^ 




We may notice that Eq/V is formally equal to the Gaussian energy density (12, 18). 
However the single-particle energy £q and the corresponding quasi-particle mass m 
entering its expression now depends on the correlated scalar density (0^), in the 
self-consistent version. 

(i^s) and {Hi) can be calculated once the RPA 2p-lp and 2p-2p GF are known (see 
the end of the previous subsection). But to calculate the total energy we also need 
the expectation value on the RPA ground-state of the one-body operators 4>q4^ and 

Ilg-nt. which are not directly given by the RPA calculation. This is the well-known 
difficulty of RPA, even in its simplest standard form, which frequently appears in 
the context of nuclear physics. In other words the calculation of the kinetic energy 
in RPA needs further manipulations. One possible way to achieve this is to use the 
so-called charging formula [11] for the calculation of the correlation energy [i.e. the 
deviation from the mean-field energy Eq) which has been historically introduced for 
the electron gas problem. Here wc will show how to adapt the charging formula 
beyond the standard RPA, namely in the r-RPA case in the context of a quantum 
field theory. 

The idea is to introduce a Hamiltonian where the coupling constant is varying be- 
tween zero and its physical value. We thus define the auxiliary Hamiltonian : 

H'{p) =Ho + pHi^t , H'{p =l) = H . (62) 



The first thing to do is to solve the RPA problem for the H'{p) Hamiltonian. For 
this purpose, one can notice that its explicit form is given by : 

H'{p) = V {^-nl .2 + ^ + ^ 1 i^n, nl + ol ,^1 0l) + p(i73 + H,) 

+ p) {\ {cP')r V : 0t + F ^ {4>y\ . (63) 



Up to constant terms, the Hamiltonian H'{p) can be rewritten as : 



+ p{H^ + Hi) (64) 

For what concerns the solution of the H'(p) RPA problem (in practice for the cal- 
culation of the commutators and double commutators entering the RPA equations) 
one has to make the following modifications with respect to the H problem : 

H ^ H'ip) 

Ol - ol = o! + l{i-p){^'u 

Hs + H^ ^ piHs + H^) (65) 

The single-particle energy occurring in the self-consistent RPA GF will be thus 
modified according to : 

^2 2 //^2 I ^ /i2\ 



= + (66) 

where (^^)p is the scalar density in the correlated RPA ground state of H'{p). Again 

the notation (<^^)_Rp is employed : in r-RPA, it coincides with the scalar density 
calculated on the self-consistent ground-state of the H'{p) Hamiltonian, while in the 
standard RPA, it coincides with the Gaussian density, i.e. calculated on the ground 
state of Hq. The solution of the H'{p) r-RPA problem is obtained formally from 
the solution of the H r-RPA problem (eq. 46-51) by simply replacing ei by eip, 

the coupling constant b hy pb and A/i = {<Pi<Pi)r by A/ip = {(j)i<j)i)Rp calculated 
sclf-consistcntly with the H'{p) Hamiltonian. 

In the standard RPA all the expectation values of (fP' are taken on the Gaussian 
ground state. In this case the energies Sp remain identical to the Gaussian single- 
particle energies £ : 



standard RPA : {(t)^)Rp = {(f>'^)R = {(j)'^)^ , eip = ei . 

Once the r-RPA problem is solved one can calculate the RPA ground state energy 
relative to the starting Hamiltonian. Since both Hq and Hint ^-re independent of 
p and since H'{p = 1) coincides with the original H one can apply the charging 
formula. The RPA ground state energy can be obtained as : 

Erpa = Eo+ ^ {pHint)p (67) 
Jo p 

where Eq is the already calculated generalized mean-field energy. Using Wick the- 
orem with respect to the vacuum of quasi-particles with energies e^, the correlated 
part can be rewritten as : 

{pHint)p = {pHs)p + {p:H, :,^), -V^ ({cl>%^ - {<l>%y 

-V^ {{4>')r - {^%,) E (0i0l)p - {^%^ ■ (68) 

In this formula {4)'^)r is as before the self-consistent scalar density of the original H. 
{^'^)e is the scalar density on the generalized mean field vacuum (vacuum of quasi- 
particles with energy e^) in the H problem and {(t>^)ep corresponds to the equivalent 
quantity for the H'{p) Hamiltonian. The remaining expectation values noted ( .. )p 
have to be taken on the r-RPA ground-state of H'(p). The calculation of these 
latter expectation values are made using eq. (46-51) where all the quantities are now 
relative to the H'(p) problem as explained before. 

In the particular case of the standard RPA the extra term in the expression of the 
correlation energy disappears since, following eq.(57), one has : 

{pHint)p = {pH^)p + {p:Hi ■.,)p . (69) 

As we will see explicitly in the next subsection the expectation value of the normal 

ordered Hamiltonian will involve an integration over the calculated Green's func- 
tions. In the case of the standard RPA the p integration can be done analytically. 
This is not the case in the r-RPA since these Green's functions will involve the e^'s 
and the self-consistent densities Mp which depend explicitly on p. 

4.3 Results in standeird RPA in 1+1 dimension 

Single particle mode 

The RPA single particle mode up with momentum P is obtained as the solution of 
the equation see e.g. eq.(46) : 



uj% = £% + T.{E = ujp, P) 



(70) 



In the standard RPA, the densities are simply taken as A/i = l/2ei where the 
ei's are the Gaussian single-particle energies. In that case I(E,P) (cq. 48) and 
consequently P) (eq. 47) are explicitly covariant in the sense that they depend 
only on E'^ — and not on E and P separately. After a simple boost-like change 
of variables one can show that : 



I{E,P) = I{E^ - P^) 
dt 1 



/ 



{2TrY Et - P^ - Aef + 



irj 



(71) 



Consequently the RPA mode has a dispersion relation which is ujp = + P^. The 
mass M of the single-particle RPA mode is thus the solution of the equation : 

= E(£;2 _ p2 = M^) . (72) 

In one spatial dimension the RPA mass operator is finite and there is no need of 
further coupling constant rcnormalization. In figure 3 the result of the calculation 
for M in one spatial dimension is shown for various values of the dimensionless 
coupling constant as a function of s. It is apparent that for p larger than a certain 
value the RPA equation may have an imaginary solution. Such a feature, which can 
appear in RPA, simply means that the HFB ground state is unstable. Hence for that 
particular theory one has to go to a superior version of the HFB-RPA approach. 

Correlation energy 

Although the standard RPA leads, in this particular theory, to an instability in a 
certain range of coupling constants, it is however interesting to look at the expression 
of the correlation energy at least at a formal level. We divide the correlation energy 
density in three pieces 

P f{3) p(4c) p(4nc) 

J^corr J-'corr J-'corr J-'corr f'70\ 

According to eqs.(50, 67), the piece corresponding to the three-body Hamiltonian 
writes : 



E^ ^ r'dp r dP r idE pH'^s'' l{E,P) 

V Jo p J {2ttY J (27r) ^ 6 1 - ^ I{E, P) ^ 



E' -P' -m' - 7' ' ^ . (74) 

2 1 - ^I{E,P)) 



Notice that, as mentioned before, the p integration can be performed analytically. 
It is convenient to transform the energy integration into an integration on the imag- 
inary axis (Wick rotation) by making the change of variable E = iz. In the usual 
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Figure 3: Squared standard RPA mass (in unit of /x^) as a function of s. 



RPA, the loop I{E, P) and thus the whole integrand /is actually depends only on 
E'^ — P'^. After the Wick rotation it depends only onS = z^+P"^ and the momentum 
integration can be done analytically in 1+1 dimension : 

For what concerns the four-body interacting piece we start from the explicit form 
of its normal ordering with respect to the e basis : 

^' '''^ " 24F JrT72iI '^^+2+3+4 ((^1 4 4 4 + ^-1 b-2 5-3 5-4 

+ 4 h\ h\ h\ 6-4 + 4 6| 6_2 5-3 5-4 + 6 5l 5| 5_3 5_4)^ (76) 

where summation over repeated indices is now omitted. As explained in subsection 
4.1, these matrix elements can be evaluated as an energy integral of two-particle 
Green's functions whose explicit expressions are given in the appendix. Noticing 
that all the matrix elements (5^5l^5^5) are identically zero, it is convenient, after 
standard manipulations, to split (: i?4 ig) into two pieces : 

(: iJ4 = (: ■■e)^"'^ + {: ^ :.)("^^ 



-2,5i+3<52+4(l + 2{b\bi)) 



b f idE ^^Er,^ 



J W) 



2AV J (27r) 



{: H4 ■.e)^'"''^ = 4(b\ bl 6-36-4) 



6 /■ idE ^^E^-i 



r idE 



2AV J (27r) 



(77) 



where the suffices (c) and (nc) stand for covariant and non covariant in a sense to 
be discussed just below. It is important to notice that the above result remains 
valid even in the case of the r-RPA where the occupation numbers (6^6) do not van- 
ish. Using the results of the appendix and the charging formula, the corresponding 
contributions to the correlation energy can now be obtained : 



where the p integration can be again performed analytically, in standard RPA. 
F{E , P , p) has the explicit expression : 



F(E,P,p) = F{E'-P\p) = —[ '^^ —+ , ^ „ 

^ ^ ^ ' - P^I(E,P) {l-^I{E,P)f 

[e^ -P^ - - '-— ' ^ . (80) 

V 2 1- £^I{E,P)) ) ^ ' 

The contribution £^iorr is explicitly covariant in the sense that the integrand depends 
only on E"^ — P^ and the trick of eq.(75) can be applied again. For what concerns 
E^rr \ the explicit calculation is more delicate since 

r dk^ dko ^^'^^ (P-ki- h] 1 

7«(S,P)= / ^ ^ , (81) 

^ ' ' J (27r)<^ 2£i2£2 E-ei+e2 + ir]' ^ ' 

appearing in eq.(79), depends separately on E"^ and P^ . However, as it is familiar in 
nuclear physics, Ectrr^ which involves expectation values of the b^b^bb^s is of higher 
order in the Y amplitudes than Ectrr involving 6666 ground state matrix elements. 



Figure 4: Diagrammatic view of the correlation energy. In the actual calculation, 
the RPA one-particle propagator (thick line) is replaced by the mean-field one. 



Indeed, it can be checked analytically that to leading order in the interaction (i.e. 
replacing F{E,P,p) by a constant value 6^/24), E^f* identically vanishes. 

Hence, we find that the correlation energy contains a piece, E^rr\ which is mani- 
festly non covariant even in the standard RPA . This problem has not been pointed 
out before, since, to our knowledge, the RPA correlation energy has never been 
calculated in a case of a relativistic theory for bosons. We have neglected this con- 
tribution in our preliminary numerical estimate for the reasons given just above. 
Nevertheless we give the explicit result, involving a four-dimensional integration, of 
this non covariant contribution. After some manipulations and change of integration 
variables, one obtains : 

E'^^f^ f^dp dS W2 



V 

with : 



I/jf^^l' ''is^o.^>AS.e)-lH-S))F(-S) (82) 

j^S,e) = - [ ^ , ^ ^ \ „ (83) 

•I £,^4e?+5sin20<5 + 44 

to be compared with : 

As mentioned above there is an instability of the HFB ground state against RPA 
fluctuations, which makes the correlation energy divergent. However, to have a first 
idea of the influence of the RPA fluctuations we replace in the above expressions for 
the correlation energy the RPA one-particle propagator by the mean field one. This 
is illustrated in figure 4 where a diagrammatic interpretation of the RPA correlation 
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Figure 5: Left panel : Effective potential for various values of the dimensionless cou- 
pling constant p in standard RPA. Right panel ; value of the condensate minimizing 
the standard RPA effective potential as a function of p. 

energy is shown. Adding this correlation energy to the mean- field energy Eq one 
obtains the RPA effective potential i.e. the RPA energy versus the condensate s for 
various values of the dimensionless coupling constant p. One gets a second-order 
phase transition with a critical coupling Pc = 1-8 (see figure 5). This has to be 
compared with the lattice result [9] and cluster expansion technique [7] showing a 
second order transition respectively at Pc = 2.55 and Pc = 2.45. It is fair to mention 
that the neglected non covariant contribution £^corr repulsive and will likely push 
the critical coupling constant to a higher value closer to the lattice result. Although 
this result is encouraging, it is obviously needed to go beyond the standard RPA to 
eliminate the unphysical instability mentioned above and seen in fig. 3. 

4.4 Results in renormalized RPA in 1+1 dimension 

Single particle mode 

In view of the calculation of the correlation energy we have now to solve the RPA 
problem for any value of p. The main problem is thus to determine the self-consistent 
M^^p = {4'q4''^^p- Oii6 very usual possibility is to calculate it self-consistently accord- 
ing to : 

J^Pp = j'-^^"'''"Gp{E,P) (85) 



where Gp{E,P) is the one-particle propagator for the H'{p) problem. One serious 
difficulty is that covariance is now lost in the sense that the loop integral Ip{E,P) 
and consequently the mass operator Tip{E,P) depends separately on E and P due 
to the presence of the density M in its expression. This is certainly a weakness of the 
present approach. However, as discussed above, even standard RPA seems to have 
problems with covariance so we think that this additional difficulty simply reflects 
the fact that there is a general problem of RPA with respect to covariance. Further 
work is needed to clarify this point. On the other hand, one natural possibility to 
recover covariance consists in imposing that the correct Ip{E, P) is obtained through 
its CM expression according to : 



Ip{E,P) = Ip{E^ - 



_ r dt 

~ J 2^ 



= ■ . (86) 



Lets us call f^pp = y-^l + the RPA single particle mode which is solution of 
the equation : 

nip = elp + E,(£;2 - p2 ^ Ml) . (87) 
In the quasi-particle approximation, the solution for Mpp is : 



Arp,-(^,-4), = ^. (88) 

The self-consistent equation for the density thus becomes an equation for the mass 
of the RPA single-particle mode which explicitly writes : 

Ml = ml + ^ [t^) " = "^'^ ^^^^ 
with the generalized mean-field single-particle mass given by : 

^i-.^. . ^ ((,^) - /_: . f ((.^ - (.^)) (.0) 

which follows directly from (66) . Hence we see that the equations for the generalized 
mean-field energy and for the single particle RPA mode are now coupled, due to the 
presence of the self-consistent scalar density {(fP')- They finally reduce to determine 
Mp. The procedure to solve the resulting equation at a given value of h and s is the 
following. We first solve the equation for p = 1, which gives the RPA mode mass 
M, the densities A/p, (0^) = "^^pNp/V and the mean-field mass m. Once this is 
done we solve for Mp which allows to obtain rup and £pp = ^P'^ + m^. 



We show on figure 6 the results of the calculation for the RPA mass M. We see 
that the instabihty problem has now disappeared. This is a first important success 
of the rcnormahzed RPA. 

Correlation energy 

The calculation of the correlation energy can now be done by assuming again co- 
variance in the sense explained just above. The results of eqs.(74, 75, 78, 79, 82) 
can be applied by just making in the final expressions, the replacements : 

dt 2Nu 

Hp 

dt 1 2Nu 



j{s,9) jp{s,e) = - [ ^ , ^ 



(92) 



In figure 7, we show the effective potential for various values of the dimensionless cou- 
pling constant p, again neglecting the non covariant piece (82) for reasons explained 
above. For p below ~ 2 there is only one minimum at s = i.e. corresponding 
to a symmetry unbroken phase. Beyond this value, a weakly pronounced minimum 
starts to develop at finite s. The symmetry broken phase becomes stable at Pc — 2.3 
indicating a very weak first order transition. It is satisfying to see that the value of 
Pc has moved in the right direction towards the value given by cluster expansion [7] 
and lattice calculation [9]. It remains to calculate the non covariant contribution to 
see if it is able to transform this weak first order transition into a genuine second 
order one. Work in this direction is now in progress. 



4.5 Towards full renormalized RPA 

In the previous section we have presented a version of r-RPA which incorporates 
RPA correlations in the scalar density which are induced by the presence of a non- 
vanishing condensate. However in the usual non symmetry broken case this version 
is still equivalent to the standard RPA. To go beyond standard RPA is, as already 
stated, a difficult problem. One possibility is to introduce the dynamical mass 
operator modifying the single particle propagator. This mass operator is of the 
form : 

E(f) (t, t') = 2e^{-i T{[H, 6„] (t) , [H, bU {t')))irr 

= 4 ^ ^ ya,-l,-2,-3(-^n'^102</)3(^), </>'lV;V'i(t')))K«,-l',-2',-3' ■ (93) 

123 1'2'3' 

Applying perturbation theory on top of the already calculated one particle propa- 
gator, one obtains for the full propagator : 



G^^^(i,iO = Go,{t,t') + j dtidt2Gait,ti)J:^^\ti,t2)Ga{t2,t') . (94) 



Using a factorization approximation, one obtains : 

= ^ E E ^",-1,-2,-3 X 
123 1'2'3' 

G,{t,t') G^^^^^^t,<^t/t,t') K.,-1',-2',-3' ■ (95) 
The density Na can be in principle calculated as : 

= {^U.) = ^ lim Gr\t, t')^ f'-^ e^^"" ^'-^(E) . (96) 

t'^t+ J ZTT 

The solution of this problem, i.e. to find self-consistently the scalar densities, is both 
formally and numerically very involved. However, as shown in a separate publication 
[24] at least in the symmetry involved region of the anharmonic oscillator, this 
procedure reproduces to leading order in the correct occupation number from 
the exact SCRPA ground state wave function. 

5 Further remarks and discussions on the formahsm. 

It is well known in non-relativistic many body theory that standard RPA corre- 
sponds to a bosonization cither of pairs of fermion operators or pairs of boson oper- 
ators. As a matter of fact, the RPA has first been invented for Fermi systems and 
the bosonization of boson pair operators has appeared much later in the literature 
[19, 20]. However, in any case the bosonization of boson pair operators (which are 
themselves NOT ideal bosons) goes along the same lines as in the fermion case [21] 
and also the generalization to relativistic field theory [6] presents no particular ob- 
stacle (as a specific example, the bosonization technics has extensively been studied 
in the case of the NJL model, see e.g. [22]). In doing so, one would naturally replace 
the boson pair operators present in the definition of the Green's function (29) by 
ideal boson operators, that is for example bq^b'^^ Bqqi^ , with Bqqi^ being an 
ideal boson operator. With such a bosonization scheme the equation of motion for 
the Green's function corresponding to 29 would lead to an inhomogeneity Mah in 
the Dyson equation of 32 which contains only unit operators on the diagonal. In- 
spection of 32, however, reveals that Nah contains in addition one body expectation 
values. This fact stems from the particular feature that in our approach we do not 
bosonise but stay with the original boson pair operators. Contrary to the usual two 
body Green's function technique, the boson pair operators are, however, taken to 
be at equal times and not at two different times. This means that we single out 
the S-channel. It is not widely known that an exact integral equation for such two 
body-two time Green's functions can be established analogous to the Dyson equa- 
tion in the one-body case. The mass operator of this 'Two Body Dyson Equation' 
presents, as in the one body case, an instantaneous mean field type part and a truly 
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Figure 7: Effective potential for various values of p in r-RPA. 



dynamical part depending on the energy. Neglecting the latter leads precisely to the 
selfconsistent equation 32. It has turned out that in this way, working with pairs 
of fermion or boson operators without bosonization, allows to better respect the 
Pauli principle and in the case of fermions this has given excellent results in a se- 
ries of applications to non-trivial models where comparison with exact solutions was 
available [14, 19, 23]. There is no reason not to believe that this advantage should 
not carry over to the case of relativistic interacting boson fields. In the relativistic 
case, however, our ansatz 29 for the Green's function has the disadvantage that it 
is not manifestly covariant and this is certainly at the origin of some difficulties we 
encountered in this respect and which wc have discussed in the main text. However, 
as we have shown in section 4, even standard RPA has, in our model at least, some 
difficulties with co-variance. Further studies are necessary to elucidate this point 
and eventually to cure it. We believe nevertheless that the results we have found in 
the present study are encouraging and that we will be able to apply our theory in 
future work as successfully to relativistic field theory as we did already in the past 
for the non-relativistic many body problem. 

6 Conclusion 

In this work we have tried to make the first step in elaborating an extension of RPA 
theory, which has been very successful in the context of non-relativistic many body 
problem [13, 15, 16, 17, 23], to relativistic field theory. Applications of standard 
RPA to relativistic field theory has emerged in the recent past and proven its great 
potential interest [4, 5, 6]. The main quality of the RPA approach is to sum a 
certain class of diagrams (the rings) to all order and by the same token to restore 
spontaneously broken symmetries and to fulfill the conservation laws. The drawback 
of standard RPA is to ignore, at a certain step of its derivation, correlations in 
the vacuum (the quasi-boson approximation). This often entails a rather strong 
overbinding of the ground state. To avoid this approximation is the aim of the afore 
mentioned extension of RPA leading to the so-called Self- Consistent RPA (SCRPA). 
An intermediate but considerably less complicated version of this theory is the so- 
called renormalized RPA (r-RPA) where, with respect to standard RPA, only the 
occupation numbers are modified due to ground state correlations. It is this latter 
version which we have tried to develop here in the context of relativistic field theory 
with application to the cp'^ theory in 1 + 1 dimension. We have studied the transition 
to a symmetry broken phase in varying the coupling constant. We have found a very 
slight first order phase transition and concluded that it will turn to second order, as 
it is expected in this model, once further correlations of the SCRPA are included. 
This opinion stems from the fact that going from the mean field theory (Gaussian 
approximation) to the r-RPA solution the first order character of the transition 
has been very much attenuated. We also point out a certain number of difficulties 



with the extension of RPA to relativistic field theory for bosons. This concerns 
for instance the fact that the approach is not manifestly covariant. Although the 
standard RPA yields at the end a covariant solution for the single particle mode, 
surprisingly we found that it has difficulties for the calculation of the correlation 
energy with respect to covariance. Apparently this had not been noticed before. At 
the r-RPA level we find that covariance is violated already for the single-particle 
mode and we have to restore it by an ad-hoc but natural assumption. It will be an 
interesting further study whether SCRPA inherently violates covariance or whether 
this is due to the approximations we have been forced to introduce. Another open 
question to be studied in the future concerns renormalization. In the present model 
study this difficulty was absent since the ^p^ theory in 1 -|- 1 dimensions is super 
renormalizable. However, in the general case, this problem has obviously to be 
mastered. Finally a detailed comparison of the diagrammatic content of the RPA 
and the cluster expansion [7], in the context of relativistic bosonic theories with a 
broken symmetry, would be certainly of great interest. 

In short we have applied for the first time an extension of RPA theory, which turns 
out to be successful in the non relativistic many body problem, to a relativistic but 

schematic field theoretic model. Although some problems are still present, we believe 
that our results are quite encouraging. Studies for the resolution of the remaining 
problems are under way. 
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8 Appendix 

In this appendix we list the explicit expressions for the Green's functions. For this 
purpose we introduce various quantities : 



1 1 EfsAfp + ep'Np' 



2e^ 2£qi E — — EjS' + ir] 



1 1 epNn + e'pMp' 



(97) 



We also introduce the loop integrals 



V 



UE) = l(^){E)+li'){E)+li'){E). (98) 
In particular for a corresponding to the momentum P one has the explicit expression: 

'ei + £2 El Ml + 82^2 



dki dk2 
£1 - £2 eiM - £2 A/2 



2 £162 E^ - (£1 + £2)^?;?? 



2£i£2 ^2 _ (^^^ _ £2)^ 

For the one-particle Green's functions one obtains : 



(99) 



E +ea + 



2£a 

-E+e^ + 5^ 



2£a 



^Go,{E) 



Ga-a'{E) — 5a,a' "} 2 ^ Ga{E), 



4£2 



G_Qt-a'(-^) 
G_Q,ta't (-E) 

where the full propagator is : 

The mass operator being given by : 

b'^s^ Ia{E) 



2 1-IUE) 



(100) 



(101) 



(102) 



For what concerns the 2p-lh and 2p-2p Green's functions we introduces indices i to 
label the destruction (creation) operators: 1 = /3, /3'(/3^, /3'^), 2 = (/?, —P'^)sym{{P^ , —P')sym) 
and 3 = -^t^ -/3'^(-/3, -/?')• The results are : 
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